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Discharge current, terminal voltage, and mass bit measurements are performed on a 
coaxial gallium electromagnetic thruster at discharge currents in the range of 7-23 kA. It 
is found that the mass bit varies quadratically with the discharge current which yields a 
constant exhaust velocity of 20 km/s. Increasing the electrode radius ratio of the thruster 
from to 2.6 to 3.4 increases the thruster efficiency from 21% to 30%. When operating with 
a central gallium anode, macroparticles are ejected at all energy levels tested. A central 
gallium cathode ejects macroparticles when the current density exceeds 3.7 xlO 8 A/m 1 2 . A 
spatially and temporally broad spectroscopic survey in the 220-520 nm range is used to 
determine which species are present in the plasma. The spectra show that neutral, singly, 
and doubly ionized gallium species are present in the discharge, as well as annular electrode 
species at higher energy levels. Axial Langmuir triple probe measurements yield electron 
temperatures in the range of 0.8-3. 8 eV and electron densities in the range of 8 x 10 20 to 
1.6 x 10 21 m -3 . Triple probe measurements suggest an exhaust plume with a divergence 
angle of 9°, and a completely doubly ionized plasma at the ablating thruster cathode. 


Nomenclature 


A = vapor pressure constant; area, m 2 
B = curve fit coefficient; vapor pressure constant 
I = current, A 

I em = electromagnetic impulse bit, N-s 

k = Boltzmann constant, J/K 

L = inductance, H 

m = ion mass, kg; mass bit, kg 

P = power, W 

P v = vapor pressure, Pa 

R = radius, m 

T = thrust, N 

T e = electron temperature, K 

u e = exhaust velocity, m/s 


Vorc = arc voltage, V 

V p = plasma potential, V 

Z = charge state; impedance SI 

a = curve fit coefficient; charge parameter 

j3 = ablation exponent 

T = particle mass flux, m -2 s -1 

e 0 = permittivity of free space, F/m 

r] = thruster efficiency 

\d = Debye length, m 

fi Q = vacuum permeability H/m 

X = non-dimensional potential 

ip = integral of current action, A 2 -s 
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I. Introduction 


T he gallium electromagnetic (GEM) thruster 1-3 is an electromagnetic thruster that ablates central elec- 
trode material for propellant mass. Gallium appears to offer numerous advantages over propellants 
currently in use. The non-toxicity and high liquid temperature range (30 °C - 2204 °C) provide numerous 
system advantages, while the high atomic mass and low ionization potential (5.99 eV) predict low frozen 
flow losses. By ablating the gallium electrode to provide discharge propellant, the electrode ablation and 
sheath losses inherent in MPD thrusters serve the useful function of vaporizing the propellant rather than 
simply introducing more heat that must be disposed of by a cooling system. In addition, mass injection is 
controlled by the arc, eliminating the need for pulsed gas valves. 

In this paper we present voltage, current, and mass bit measurements taken at discharge current levels 
in the range of 7-23 kA. Emission spectroscopy measurements are undertaken to examine the atomic and 
ionic species present in the discharge. Axial triple probe measurements are taken to determine the electron 
temperature and density along the center-line of the the thruster. 


I. A. Basic Scaling Relations 

The thrust of self-field electromagnetic thrusters follows a quadratic dependence 4 on the discharge current 


T= -L'l 2 , where V = 2 
2 


£ fin— + j 
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The exhaust velocity is found by dividing the thrust by the mass flow rate 


T 

u e — . 
m 


(i) 


( 2 ) 


For ablative- type thrusters, it is convenient to write the mass flow rate as following a power law dependence 
with the discharge current 


m oc I& (3) 

where (3 is an ablation exponent. Note that for gas-fed thrusters, the mass flow rate is independent of the 
discharge current and (3 = 0. Substituting Eqns. (1) and (3) into (2), we find 

u e oc J 2- * 3 (4) 

The power across the discharge is given by 

P = IV arc oc imu 2 (5) 

Prom which we find 

V arc oc I 3 ^ (6) 

Therefore, by examining the V-I characteristics and mass bit as a function of current, we can obtain a 
general understanding of thruster performance. The performance of the Teflon PPT, 5 Teflon MPD, 6 and the 
graphite MPD 7 has been investigated. For all three thruster types the discharge current, arc voltage, and 
thrust were measured as a function of input energy. The exhaust velocity was found by dividing the impulse 
by the measured mass bit. In both Teflon- fed thrusters, the propellant is ablated when an arc is struck 
between coaxial metal electrodes. The graphite MPD is similar to the GEM thruster in that the electrode 
material is directly ablated as the propellant. 

In order to compare the performance of ablative thrusters with varying pulse lengths and discharge 
current levels, the integral of the current action is introduced 

+ =r>(t)dt (?) 

JO 

The V-I characteristics and ablated mass as a function of ip are shown in Fig. 1. Both plots yield linear 
results for the graphite and Teflon MPD. Prom both equations (3) and (6), this yields (3 = 2, or a constant 
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exhaust velocity. Indeed, it was found that the exhaust velocity was invariant with the discharge current, 
with values of 28 km and 14 km/s, for the graphite and Teflon MPD, respectively. The PPT has an ablation 
exponent value of (3 = 1.5, and it was found that the exhaust velocity increased with discharge current. 5 
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Figure 1. Voltage-current and mass bit vs. ijj plots for previously developed ablative electromagnetic thrusters. The 
ablation exponent is /3 = 1 for both the graphite and Teflon MPD, while for the Teflon PPT (3 = 1.5. 


II. Experimental Apparatus 


II. A. Testing Facilities 

Tests were conducted at the University of Illinois Electric Propulsion Laboratory on the coaxial GEM thruster 
shown in Fig. 2. The stainless steel annular electrode has outer and inner diameters of 57.2 mm and 22.6 
mm, respectively. A 8.7 mm diameter hole drilled through the center of the boron nitride insulator allows the 
gallium to make electrical contact with a brass standoff, which in turn is connected to a high-current pulse 
forming network described below. The thruster was mounted into a 5-way vacuum cross, which is attached to 
the main vacuum chamber. A baffle was inserted into the 5-way cross to condense and capture ejected gallium 
particles before they reach the main chamber. Vacuum was maintained by a TP HI 500 turbomolecular pump, 
which is backed by two Roots blowers and two Kinney mechanical displacement pumps. The base pressure 
of the facility during testing was 5 x 10 -5 torr. 



Figure 2. Coaxial GEM thruster electrodes, and boron nitride insulator. 


A ten-section, 5-line pulse forming network (PFN) was used to power the thruster. Each of the 5 sections 
consists of ten 88 fiF capacitors in series with 85 nH inductors. The 50 capacitors have a maximum charging 
voltage of 1.2 kV, which yields a bank energy of 3.2 kJ. A copper sheet with punched holes was located over 
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the capacitors to provide a common ground (ground plane). Brass standoffs with a height of 3.8 cm were 
machined to provide spacing between the stage inductors and the ground plane. Two 9 cm long, 0.32 cm 
(1/8) diameter copper rods yield a section inductance of 85 nH. The five PFN lines were each connected 
in parallel using 2.5 cm wide copper strip. To minimize stray inductance in the charging circuit, the PFN 
was mounted close to vacuum chamber, as shown in Fig. 3. A 7.0 cm wide strip of copper sheet formed the 
transmission line to the thruster. The capacitors were charged using a Hipotronics DC power supply. 



Figure 3. 3.2 kJ Pulse Forming Network and 1.5 m 3 vacuum tank. 


III. Experimental Results 


II I. A. Voltage and Current Measurement 


Measurements were obtained for capacitor charging voltages in the range of V 0 = 150-350 V. A Pearson 
4418 current monitor and Tektronix P6015A voltage probe were used to measure the discharge current 
and terminal voltage, respectively. Representative traces are in Fig. 4. Note that the PFN impedance is 
well-matched to the arc impedance, minimizing current reversal. 
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Figure 4. Representative voltage and current traces taken at 290 J. 


A plot of the arc voltage vs. the discharge current is shown in Fig. 5 for an anode and cathode radius 
of 15 and 4.35 mm, respectively. The voltage and current data is taken at r = 50 (i s. The arc impedance 
varied from 6-7 mfi, decreasing slightly with energy level. The arc voltage varies linearly with the current, 
resulting in an ablation exponent of /? = 2. 
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Figure 5. Plot of arc voltage as a function of discharge current. 


III.B. Mass Bit Measurement 


The mass ablated per pulse was measured by weighing the inner and annular electrode after a prescribed 
number of shots using a Mettler AG245 measurement scale with a resolution of 0.1 mg. Several hundred 
shots were fired which liberated at least 10 mg, leaving the error due to instrumentation at less than 1%. 
Tests were repeated 4 times at each energy level, and the ablated mass measurements were found to be 
repeatable. No visible oxidation was seen on the gallium surface before or after testing. The mass of the 
annular stainless steel electrode decreased slightly after each test run. However, the total amount of mass 
lost (i.e. 0.3 mg after 250 shots) was small compared to the total mass of ablated gallium. 

Measurements were initially obtained using a 5 mm central gallium anode. It was found that in this 
configuration gallium sputtered from the central electrode leaving macroparticles on the insulator and outer 
electrode (Fig. 6). Macroparticle ejection can be eliminated by reducing the current density. This was first 
attempted by lowering the discharge current, and when this failed, by increasing the diameter of the gallium 
electrode from 5, 7.1 and then 8.7 mm. This proved unsuccessful, so the polarity of the thruster was then 
reversed, and the measured mass bit dropped by an order of magnitude and no sputtered gallium particles 
were observed when using the 8.7 mm diameter gallium cathode. However, when the current exceeded 
roughly 22 kA (j = 3.7 xlO 8 A/m 2 ), gallium once again sputtered from the central electrode. This set the 
upper limit of testing at a capacitor charging voltage of V Q = 350 V (290 J). The mass bit vs. tp is shown in 
Fig. 7. The linearity of the plot is consistent with the V-I plot and also yields a mass ablation coefficient of 
0 = 2 . 
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Figure 6. 


Macroparticles sputter onto the insulator and annular electrode when utilizing a central gallium anode. 


The exhaust velocity can be estimated through use of the electromagnetic impulse bit 

Iem = f Tdt = 


(8) 
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Current Action Integral ^ ? A 3 -s 


Figure 7. Measured mass bit vs. ij> parameter. 


The exhaust velocity is found by dividing the impulse bit by the mass bit, m 

_ Iem _ 1 L'i/j 
Ue m 2 m 

The slope in Fig. 7 is proportional to exhaust velocity, which has values of 17 and 20 km/s for radius ratios 
of 2.6 and 3.4, respectively. 

III.C. Langmuir Triple Probe Measurements 

III.C.l. Probe Construction 

The triple Langmuir probe 8 has been used successfully used to measure the electron temperature and density 
in MPD thrusters. 9-11 The triple probe (shown schematically in Fig. 8) consists of three separate electrodes 
configured such that the electron temperature and density can be reduced from the measured currents I 2 
and current 73 . 



Figure 8. Triple Langmuir probe circuit. 


The triple probe used in this study consisted of three 0.381 mm (0.015") diameter tungsten wires fed 
through a four -bore alumina tube with an O.D. of 4.7 mm. The exposed wire length was 4.5 mm yielding 
a probe surface area of 5.5 zb 0.2 mm 2 . The tungsten wires are connected to RG58 coaxial cables which are 
fed to the oscilloscope located in a shielded screen room. The applied voltages of Vd 2 = 5 V and Vaz = 10 V 
were supplied using a DC power supply floated through use of an isolation transformer. 
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III.C.2. Probe Theory 

Use of a triple probe assumes (a) a Maxwellian electron energy distribution, (b) a thin ( r p » A d) , 
collisionless » A d) ion sheath, and (c) an ion sheath thickness smaller than the separation between 

probes ( d s » Ad)- The ion sheath is on the order of the Debye length, Ad, which is given for a quasi- neutral 
plasma by 


A d 



( 10 ) 


Assuming T e = 2 eV and n e = 10 21 m -3 , yields a Debye length of 0.3 microns. Therefore 1 >> r v j\n and 
the thin sheath requirement is met. The ratio of the ion-ion and ion-electron mean free paths to the Debye 
length are 150 and 24000, so the sheath can be assumed to be collision-less. The separation of the probes is 
approximately 0.5 mm, which is much greater than the Debye length. 

The current flowing into the probes may be written as 


h = AJ e exp(-xi) - A Ji(B + XiT 


(ii) 


h = AJ e exp(— * 2 ) - AJi(B + X 2 ) c 


I 3 = AJ e exp(— x 3 ) - AJi(B + X3)“ 

where the dimensionless potential x for each probe is given in terms of the plasma potential V p . 

e\V-V p \ 

X kT e 

The electron and ion saturation current current densities are given by 

_ ( kT R N 1/2 

J p — Tip e 


J * — e 


\ 2i xm e 
kT e ^ 1/2 


‘hvnfii 


The coefficients a and B are found through use of the Peterson/ Talbot curve 12 

2.9 


a 


B = 1.5 + 


ln(r p /AD) + 2.3 

Ti 


rp \ 3/4 

+ 0.07 ( -^r 1 - 0-34 


ZjT P 


0.85 + 0.135 ( In ( -p-)) 


Equations (10)- (18) result in a non-linear system of equations that yield T e , n e , and V p . 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


III.C.S. Triple Probe Results 

Measurements were taken on the center line of the thruster downstream of the cathode tip, at 2, 6, 12, and 
18 cm. All measurements were taken at a discharge current of 23 kA. Measurements were averaged over 
several shots and the probes were manually cleaned in between experiments. The measured currents I 2 and 
current Is are shown in Fig. 9 for applied voltages of V ^2 = 5U and Vds = 10V\ The applied voltages where 
measured during each pulse and found to vary by less than 10% of the nominal value during the steady 
portion of the pulse. 

The axial variation of the peak electron temperature and density are shown in Figs. 10 and 11, with 
values at an axial location 2 cm downstream of the inner electrode of 3.8 eV and 1.6 xlO 21 m -3 . The 
sources of experimental error include: uncertainties in probe surface area and probe end effects (10%), probe 
misalignment (5 %), the presence of doubly ionized species (20%). The error for the electron temperature is 
thus estimated to be 15% with a 35% error on the number density. 
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Figure 9. Collected triple probe currents taken 2 cm from cathode tip. 
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Figure 10. Center line electron temperature at four axial probe positions. 



Figure 11. Center line electron density at four axial probe positions. 
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Extrapolating the temperature data to the face of the cathode (z = 0 cm) and assuming an expansion 
ratio of 100 yields an electron temperature of 5 eV and an electron density of 10 23 m~ 3 . At this temperature 
and density range, the plasma is completely doubly ionized. This results in a Larmor radius of 15 mm, 
which is approximately the radius of the annular electrode This is suggests that the bulk of the ions will 
be able to drift out of the discharge region without being lost to the outer electrode. An estimate of the 
plume divergence can be made by assuming that the product of the particle density, exhaust velocity, and 
plume area is a constant ( n e u e A = const). The diameter of the expansion fan is taken to be 5 cm at a 
downstream distance from the cathode surface of two centimeters (Fig. 12). This yields a divergence angle 
of tan- 1 (2.5/16) =9°. 


ijg =20 km/s = const 


i 



Figure 12. Plume divergence for GEM thruster. 


III.D. Emission Spectroscopy 

We acquired spectroscopic data from the discharge to determine which species were present. These spectra 
were obtained using a 25 mm focal length spectrometer with 5 micron inlet slit, yielding 0.7 nm resolution 
over the range of 220 nm to 520 nm. Emission was collected using a 12 mm diameter f/2 quartz piano 
convex lens which was focused onto the end of a 400 micron diameter silica fiber. The spectrum was time- 
integrated over several discharge pulses and spatially integrated over a broad region. The output of the 
spectrometer was wavelength-calibrated using a mercury lamp and corrected for relative spectral efficiency 
using a deuterium lamp. 

The emission spectra were obtained at energy levels ranging from 85-290 J. Emission from neutral, singly, 
and doubly ionized Ga species were observed in the discharge. Iron, removed from the outer electrode, was 
also seen in the spectra obtained. Spectra from the 85 and 290 J cases are shown in Figs. 13 and 14. They 
are largely similar, however the 290 J case exhibits several iron lines as well as additional Ga III lines. A 
full listing of the observed gallium species is presented in the Appendix. 



Figure 13. Emission spectra taken at energy level of 85 J. 
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Figure 14. Emission spectra taken at energy level of 290 J. 


IV. Analysis 


IV. A. Efficiency Calculation 

The efficiency of electric thrusters is defined as the ratio of the jet kinetic power to the input power 

raui 

V= Wn 

An alternative form can be obtained by utilizing the dynamic impedance 


(19) 


Zdyn — ^ u e (20) 

and by substituting P = IV arc in equation 19. One obtains 

(21 , 

Z arc 8rnZ arc 

where Z arc = Zd yn + An estimate of the thruster efficiency can be made from the discharge current, 

arc voltage, and mass bit. Tests were initially conducted on a thruster possessing a radius ratio of 2.6, 
which yielded a computed efficiency of 21%. Increasing the diameter of the annular electrode resulted in an 
increase in the efficiency to 30%. Assuming the electrical characteristics and mass bit are roughly constant 
for a given discharge energy, the efficiency rises to 40% at a radius ratio of 5 (Fig. 15) with an exhaust 
velocity of 23 km/s. 



Figure 15. Exhaust velocity and thruster efficiency as a function of electrode radius ratio. 
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IV. B. Scaling Relations 

A review of MPD performance has been presented in Ref. [2]. Two parameters of particular interest are the 
ion Larmor radius Rl , and the charge carrier parameter a, which represents the fraction of the discharge 
current that can be carried by the number of ions present 


Ze dm 
ml dt 


An equivalent Larmor radius can be defined as: 


( 22 ) 


Rl = 2ttR— 

Zeu e ji 0 l 

where fi 0 is the magnetic permeability of free space and R is the average radius of the coaxial thruster. The 
Larmor radius and charge carrier parameter are tabulated for different propellants. It is assumed that both 
the ions and electrons carry the current, so the charge carrier parameter is multiplied by 2 in the Table IV.B. 
Additionally, for both ablative thrusters (the GEM and Carbon MPD) the charge state is taken to be Z = 2. 


Table 1. MPD Thruster performance for various propellants 


Propellant 

At. Wt. 

Current, kA 

dm/dt, gm/s 

Isp, s 

Eff. 

Charge Parameter, a 

Larmor Rad., cm 

Hydrogen 

1 

20 

0.5 

16000 

0.52 

4.8 

1.3 

Carbon 

12 

24 

0.6 

2780 

0.07 

0.7 

0.4 

Nitrogen 

18 

18 

4.0 

4000 

0.37 

3.2 

5.0 

Argon 

40 

20 

6.0 

2800 

0.20 

1.4 

8.8 

Gallium 

70 

23 

4.4 

2000 

0.30 

1.4 

1.5 

Xenon 

131 

20 

6.0 

1400 

0.10 

0.4 

14.4 


The results from Table IV.B suggest that efficient operation is reached when a > 1 and the Larmor radius 
is smaller than the thruster channel. The efficiency tabulated for the GEM thruster (30%) corresponds to 
an electrode ratio of 3.4. As noted earlier, increasing the electrode ratio to 5 would yield an efficiency of 
40%, which compares favorably to the thrusters shown. 


IV. C. Mass Flow Rate Prediction 


As noted previously the mass scales quadrat ically with the discharge current. An expression for the mass 
flow rate as a function of current can be found as follows. An estimate of the evaporation rate can be made 
by using the Langmuir equation, which expresses the flux of atoms released from a surface in equilibrium in 
terms of the surface temperature and vapor pressure 13 


r = 


Py 

y/27T mkT s 


(24) 


The vapor pressure P v (in Pa) and surface temperature T s (in K) are related by the Clausius-Clapeyron 
relation 


log (P v ) = 5.006 + A - ^ (25) 

s 

where the constants have values of A = 5.52 and B = 13798, respectively for gallium. 14 The pressure due 
to the pinch effect is given by 


, _ 2 Mo/ 2 

av3 3 8 ?r 2 r2 


(26) 


where the factor of 2/3 has been added to account for the parabolic distribution of pressure over the cathode 
surface (Fig. 16). Setting the vapor pressure equal to the average pinch pressure (P v = P aV g) yields an 
expression for the mass flow rate in terms of the discharge current. A comparison of the calculated vs. 
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experimental mass flow rate is shown in Fig. 17. It can be seen that the mass flow rate predicted using the 
pinch model is consistently higher than the measured value. This may be due to a lower than expected liquid 
gallium surface temperature, due to heat conduction into the liquid. An additional effect is a non-uniform 
current density, which could reduced the 2/3 factor and also the ablated mass. 


Arc Column 1 



Figure 16. High current pinch discharge with evaporating gallium electrode. 



V. Summary 

The results obtained can be summarized as follows: 

• Experiments were conducted at discharge currents spanning a range from 7-23 kA utilizing a central 
gallium electrode. When using a central gallium anode, macroparticles were ejected at all energy 
levels tested. Macroparticles were ejected using the central gallium cathode when the current density 
exceeded 3.7 xlO 8 A/m 2 . 
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• The measured mass bit follows a quadratic dependence on the discharge current. This yields a constant 
exhaust velocity of 17 km/s and 20 km/s for radius ratios of 2.6 and 3.4, respectively, and a predicted 
velocity of 23 km/s for a radius ratio of 5.0. 

• The performance increases as the radius ratio of the thruster increases. An efficiency of 30% was 
obtained with a radius ratio of 3.4, and would be raised to 40% for a radius ratio of 5.0 

• Triple Langmuir probe measurements along the center-line of the thruster yield a peak electron tem- 
perature and density of 3.8 eV and 1.6 21 m -3 two centimeters from the cathode tip. The axial mea- 
surements suggest a well collimated exhaust plume with a divergence angle of 9°, and fully doubled 
ionized plasma at the cathode surface. 

• Emission spectroscopy was used to determine the species present in the discharge as a function of 
discharge current. Ga I-III lines were present at all energy levels. The number of doubly ionized 
gallium and iron (from the outer electrode) lines increased as the current increased. 

Appendix 


Table 2. Neutral and ionic gallium species present in 22 kA arc discharge. 


Specie 

Wavelength, nm 

Upper Energy Level, cm 1 

Lower Level 

Upper Level 

Gal 

223.592 

45071.75 

3d 10 4s 2 4p 

4s 2 7d 

Gal 

229.786 

44332.27 

3d 10 4s 2 4p 

4s 2 8s 

Gal 

233.824 

43580.44 

3d 10 4s 2 4p 

4s 2 6d 

Gal 

245.007 

40802.86 

3d 10 4s 2 4p 

4s 2 5d 

Gal 

260.734 

38341.72 

3d 10 4s 2 4p 

4s4p 2 

Gal 

294.417 

34787.85 

3d 10 4s 2 4p 

4s 2 4d 

Gal 

403.298 

24.788.53 

3d 10 4s 2 4p 

4s 2 5s 

Gal 

417.204 

24788.53 

3d 10 4s 2 4p 

4s 2 5s 

Gall 

231.727 

113842.30 

4s4p 

4s4d 

Gall 

235.532 

156327.39 

4s4d 

4s7f 

Gall 

251.403 

147485.46 

4p 2 

4s 5f 

Ga II 

277.997 

106662.37 

4s4p 

4s 5s 

Ga II 

299.265 

151923.93 

4s5p 

4s 8s 

Ga II 

337.595 

137342.57 

4p 2 

4s4f 

Ga II 

347.208 

147520.34 

4s5p 

4s6d 

Ga II 

425.593 

137332.28 

4s4d 

4s4f 

Ga III 

208.948 

209148.01 

3d 10 5p 

3d 10 7s 

Ga III 

236.866 

227637.2 

3d 10 4f 

3d 10 7g 

Ga III 

241.787 

185432.59 

3d 10 4d 

3d 10 4f 

Ga III 

242.454 

185432.59 

3d 10 4d 

3d 10 4f 

Ga III 

351.737 

189187.72 

3d 10 5p 

3d 10 5d 

Ga III 

438.184 

208254.14 

3d 10 4f 

3d 10 5g 

Ga III 

486.303 

161304.40 

3d 10 5s 

3d 10 5p 

Ga III 

499.393 

160765.56 

3d 10 5s 

3d 10 5p 
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